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We present an analytic study of domain-wall statics and dynamics in ferromagnetic nanotubes with spin-orbit- 
induced Dzyaloshinskii-Moriya interaction (DMI). Even at the level of statics, dramatic effects arise from the 
interplay of space curvature and DMI: the domains become chirally twisted leading to more compact domain 
walls. The dynamics of these chiral structures exhibits several interesting features. Under weak applied currents, 
they propagate without distortion. The dynamical response is further enriched by the application of an external 
magnetic field: the domain wall velocity becomes chirality-dependent and can be significantly increased by 
varying the DMI. These characteristics allow for enhanced control of domain wall motion in nanotubes with 
DMI, increasing their potential as information carriers in future logic and storage devices. 

PACS numbers: 75.78.Fg, 75.60.Ch, 75.70.Tj 


I. INTRODUCTION 

In recent years, ferromagnetic nanostructures featuring nar¬ 
row and stable domain walls (DWs) have been in the spotlight 
of experimental and theoretical research, with an overarching 
aim to achieve more compact spintronic logic and memory 
devices.H3 In particular, numerous efforts have been focused 
on DWs in ferromagnetic nanowires,®^ nanotubes,and 
thin films with perpendicular magnetic anisotropjl2S®featur- 
ing the Dzyaloshinskii-Moriya interaction (DMI).l5iEl Here 
we report striking effects arising from the interplay between 
space-curvature and DMI in ferromagnetic nanostructures, 
leading to narrow and stable DWs controllable, efficiently and 
reliably, by means of electric current and magnetic field. 

Curvature effects play a significant role in various fields 
of physics, and are attracting increasing attention in con¬ 
densed matter, particularly in nanomagnetism. The simplest 
system with curvature where DW dynamics can be consid¬ 
ered is a magnetic nanotube. Furthermore, thin ferromag¬ 
netic nanotubes have attracted recent attention from experi¬ 
mentalists o wing to a number of technologically advantageous 
properties, ISEEEl including enhanced DW stability under 
strong external fields, al lowing for higher DW velocities com¬ 
pared to flat geometries pSESHI] increased DW velocities under 
electric current pulses® and the possibility of swi tchin g chi¬ 
rality in vortex DWs through magnetic field pulses.®^ 

We show that in thin ferromagnetic nanotubes, the 
DMI induces qualitatively different effects to those found 
in flat nanostructures, such as thin films and rectangular 
nanowires.® In nanotubes, DMI causes the domains to be¬ 
come twisted, with helical lines of magnetization as in Fig. [T] 
forcing the DWs to become narrower. In contrast, in rectangu¬ 
lar nanowires with DMI, the magnetization far from the DW 
remains parallel to the wire axis while DWs become broader. 
This sharpening effect of DMI in nanotubes can enable sub¬ 
stantial downscaling in future nanodevices. 

We further demonstrate that in a certain thin-nanotube 
regime specified below, DWs exhibit perfectly stable motion 
under an applied electric current, propagating without any dis¬ 



FIG. 1. (Color online) Domain wall profile in a thin nanotube with 
Dzyaloshinskii-Moriya interaction. The magnetization lies tangent 
to the nanotube surface. 


tortion. The adiabatic spin-transfer torqud®® is absent in the 
spin dynamics equations, and the non-adiabatic term takes the 
form of the adiabatic one. 

Complimentarily to the current, a magnetic field along the 
nanotube triggers a rich dynamical response in the magnetiza¬ 
tion texture. We show that the DW velocity becomes strongly 
dependent on polarity and chirality,® and can be significantly 
enhanced by DMI, which is favorable for memory applica¬ 
tions. Moreover, the onset of magnonid® breakdown, imped¬ 
ing DW transport at high fields, can be efficiently suppressed 
by DMI. 


II. STATICS 


We consider a ferromagnetic nanotube with inner radius R 
and thickness w. In the thin-nanotube regime w R, the 
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micromagnetic energjEHl with DMP®II takes the form 

i?(m) = J d^r l^lVnip + X [l — (m • e^)^] 

+ DM^m • (V X m) + (m . ep)^| , (1) 

where the integral runs over the volume of the nanotube, 
A is the exchange constant, K is the easy-axis crystalline 
anisotropy, D is the DM1 constant. Mg is the saturation mag¬ 
netization, and po is the magnetic permeability of vacuum. 
The cylindrical-coordinate unit vectors e^, Op, and are 
shown in Fig. [T] 

The last term on the right-hand side of Eq. Q represents 
the shape anisotropy stemming from the thinness of the nan¬ 
otube. In nanotubes with radius R much larger than the mag¬ 
netostatic exchange length this term forces the 

magnetization to lie nearly tangent to the surface. In this case, 
the unit vector of magnetization may be described by its ori¬ 
entation Q{z, p, (j)\ t) in the (z, 0)-tangent plane; 

m = cos 0 - 1-00 sin 0 . ( 2 ) 

Substituting Eq. Q into Eq. Q and introducing the dimen¬ 
sionless coordinates s = r/R, Q = z/R, ^ = p/R, anisotropy 
K = KR^/A, and DMI constant ry = DM^R/{2A), we 
thereby obtain the expression for the dimensionless energy 


£ = E/{2AR): 

£ = J d^s{ei{Q) + e2iQ)), (3) 

where the energy densities ei and £2 are given by 

£i = i|ac0|'-y(0) + |(i + K), (4) 

£2 = I \d^Q + ■ (5) 

The “potential” V, which appears in ei, is given by 

V{9) = cos2{9 — S), (6) 


a= [(l + Kf+4p^l^ , tan2S = --^, (7) 

1 -f At 

where S is taken between —7r/2 and 7r/2. Below we shall see 
that (5 determines the orientation of the twisted domains, while 
1 /a is the DW width. 

Next we look for a magnetization profile 0 which mini¬ 
mizes the energy £. The £2 term vanishes (and is thus min¬ 
imized) by taking = —77 and 900 = 0. Then 0 is of 
the form 0 o(C) ~ ~ !)■ In th® thin-nanotube limit (and 

taking At, p = 0(1)), the £i term can be simultaneously min¬ 
imized by taking 0o(C) to satisfy the Euler-Lagrange equa¬ 
tion 0Q = —V'{9q), subject to the boundary conditions that 
00(± 00 ) correspond to maxima of V (not minima). These 
maxima, given hy 9 — S + mr, describe the orientations of 
magnetization in domains. In the case of zero DMI, ie, p = 0, 
the magnetization far from the DW center is parallel to the 
nanotube axis. However, for p 0 the magnetization profile 
becomes helical, as shown in Eig. [T] 


Domain walls correspond to boundary conditions 0o(ic«) 
describing oppositely oriented domains. There are four dis¬ 
tinct DW profiles, characterized by polarity a and chirality x 
(one is shown in Eig.[T]l, for more details see Appendix [A| Po¬ 
larity determines whether the DW is head-to-head {a = 1) or 
tail-to-tail (a = —1), while chirality determines the sense of 
rotation of m with increasing so that x = crsgn0g. The 
Euler-Lagrange equation may be solved exactly to obtain 

00 = 2xarctan -f 6 . (8) 

The DW profiles may be understood qualitatively in terms of a 
mechanical analogy - see Eig.|^ We regard 9{() as the trajec¬ 
tory of a particle moving in a potential V(9) with ( playing the 
role of time. In the static case, the DW boundary conditions 
correspond to the particle approaching consecutive maxima 
of V (located at S mod tt) as C —>■ ± 00 . At times in between, 
the particle traverses the intervening potential well (this is an 
example of a so-called instanton orbit). 



FIG. 2. (Color online) Mechanical analogy: The profile ff(() may 
be regarded as the instanton orbit of a particle 0 taking infinite time 
p to move between consecutive local maxima of V{9) — he cos 9 
(at = 1 throughout). Dashed curve: With no DMI or applied field, 
the local maxima are 9 = 0 and 0 = tt, corresponding to domains 
aligned along the nanotube axis. The instanton orbit, indicated by the 
arrows, describes a tail-to-tail DW with negative chirality. Dotted- 
dashed curve: With DMI parameter p = 1 but no applied field, the 
local maxima are shifted by 5 = — vr/8, corresponding to twisted do¬ 
mains. The instanton orbit corresponds to a head-to-head DW with 
positive chirality. Solid curve: With p = 1 and applied field he = 1, 
the values ofV—he cos 0 at consecutive maxima are no longer equal; 
a specific value j -|- v, the coefficient of linear damping in Eq. GD, 
is required to ensure that particle reaches the second maximum with¬ 
out overshooting. Hashed curve: For large enough he, a maximum 
and minimum of V — he cos 0 coalesce, and the instanton orbit is 
destroyed. 

An important characteristic of the DW is its width A given 
by 1 /a, or in physical units. 


[{K + A/R?Y -F ■ 

As is clear from this expression, DWs in thin nanotubes be¬ 
come sharper in the presence of DMI, in marked contrast to 
the case of rectangular nanowires .I^DWs also become sharper 
in nanotubes with higher curvature 1 jR. 
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III. DYNAMICS 


Under an applied current, the magnetization dynamics in a 
ferromagnet far below the Curie temperatur e is d escribed by 
the Landau-Lifshitz-Gilbert (LLG) equation®^3tlII 


9m 

1h 


TT 9m ^dm 

7 H X m + am x —- J — h pJm x 

at oz 


dm 

Ih 


, (10) 


where H = —{Ms)~"^5E/5m is the effective magnetic field, 
7 is the gyromagnetic ratio, a is the Gilbert damping constant, 
J is the current along the nanotube in units of velocity, and 
P is the nonadiabatic spin-transfer torque parameter. In the 
regime 


w, 



< R 


( 11 ) 


(as considered in the static case) and for currents J satisfying 


J < 


')A 


( 12 ) 


it can be shown that m lies nearly tangent to the nanotube, 
and Eq. ( [T0| ) reduces to (see Appendix [B] for the details of this 
calculation); 


— =-m X (m X ht) - j—. ( 13 ) 

or oQ 

Here t = 1 is the dimensionless time, ht is the tan¬ 

gential component of the dimensionless effective field h = 
H, and j = is the dimensionless current. Note 

that j is proportional to nonadiabatic torque parameter P, 
whereas the current term itself assumes the adiabatic-torque 
form. This has an important consequence, as described be¬ 
low. 

Proceeding as in Eq. (|^, we write m = cos ©-t-e^ sin 0 
with 0((^, r) = 9{C,, t) + ri{^ — 1 ) to obtain 


80 

dr 


8^0 



+ V'{0). 


(14) 


We look for traveling-wave solutions of the form 0{C,t) = 
-d^d^—VT) describing axially symmetric DWs propagating with 
velocity v. Erom Eq. ( [T4l i, the profile d satisfies 

= (j - v)^)' - U'(r?) (15) 


subject to the same boundary conditions as in the static case. 
It is easy to see that the moving profile coincides with the 
static profile 0q in Eq. with velocity v = j. In physical 
units, the DW velocity is given by 

V = pj/a . (16) 

Erom Eqs. ( [TT] ) and ( [T^ , Eq. ( [T6| ) holds for velocities in the 
regime 


Thus, under an applied current, the DW propagates without 
distortion and with velocity independent of polarity, chirality 
and DMI.EIAs the velocity approaches the magne¬ 

tization acquires a non-negligible radial component, and new 
behavior can be expected to appeal®® 

Next we study the effect of an external magnetic field on 
DW dynamics. We take the field to be uniform along the 
nanotube axis. He = H^ez- In the thin-nanotube limit, the 
applied field generates an additional term in Eq. ( [T5| ): 

■d” = {j — v)^'— V'{d) — heSind , (18) 

where he = ^f^He- The boundary conditions are modi¬ 
fied so that r?(±oo) correspond to consecutive maxima of a 
modified potential, V(-0) — he cos d. In terms of the mechan¬ 
ical analogy (Eig.[^, i?(^) again describes the trajectory of a 
particle moving from one potential maximum to another, as 
above. However, the potential difference at consecutive max¬ 
ima induced by the applied field is compensated now by the 
additional (anti)damping term (j — v)^'. 
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FIG. 3. (Color online) DW velocity v vs DMI parameter 7 for differ¬ 
ent chiralities x and polarities a (he = k — 1 and j = 0). Changing 
the sign of x and 7 leaves v unchanged, while changing the sign of 
a and 7 alters the sign of v. 


Numerical solutions of Eq. ( |T8] l show that the applied field 
causes the DW velocity to depend strongly on DMI, chirality, 
and polarity, see Eig. Eor a given field strength, the veloc¬ 
ity achieves a maximum for a nonzero value of 7 , and varies 
with DMI through 7 by a factor exceeding 2. While Eq. 
cannot be solved analytically, one can develop an expansion 
in powers of he (see Appendix]^ for the details): 


• , Vl + K + a2 TTT] 2 


(19) 


As shown in Eig. this quadratic approximation is in good 
agreement with the numerical results. In the limit of no cur¬ 
rent and DMI, j = 7 = 0, it yields v = ahej (v/Sa) in accord 
with Ref. [55] or in physical units the DW velocity due to mag¬ 
netic field reads 


V = 


^7 R 

V2a y/l + KR^/A 


( 20 ) 


V < ^-^/JlqA. 


In the limit of R ^ 00 this expression reduces to a well known 
(12) result for the velocities of transverse DWs in flat nanostrips.®! 
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FIG. 4. (Color online) DW velocity v vs DMI parameter rj for dif¬ 
ferent values of the external field he (cr = -F, x = +> '^ = 1 . and 
j — 0 throughout). The solid curves are obtained from the numerical 
solution of Eq. jl 8 [l; th e dotted curves are given by the approximate 
analytical formula ip^. 



FIG. 5. (Color online) DW velocity v vs. applied field he for differ¬ 
ent values of the DMI parameter 77 and chirality x (f = -F, ft = 1, 
and j — 0). For 77 = 0, u is given by the exact linear relation 
V = ahe/a. For 77 7 ^ 0, 77 is obtained by solving Eq. numeri¬ 
cally. Curves are computed up to the critical field he. 


The dependence of the DW velocity on applied field and 
chirality is shown in Fig. It follows from Eqs. ([T9) and 
(0 that for small fields, the DW velocity is suppressed by the 
DMI. However, for larger fields and chirality x = o' sgn{rihe), 
the velocity may be enhanced by DMI (for the opposite chi¬ 
rality, V is always reduced). 

At a certain critical applied field he, a bifurcation occurs, 
beyond which the DW velocity is suppressed. In terms of the 
mechanical analogy of Fig. as he approaches he, a maxi¬ 
mum and minimum of the potential V — he cos d coalesce, the 
instanton orbit is destroyed, and the character of the traveling 
DW changes. This phenome non h as been discussed in terms 
of the spin-Cherenkov effect,l2S^ and more recently in terms 
of pulled wavefronts of the KPP equation.l^ It is straightfor¬ 
ward to obtain an analytic expression for he in terms of rj, 
shown in Fig.|^ For 77 ^ 1 + k the leading-order behavior is 
given by /ic = 1 + K — 3{cri)^^^, with c = v/l + k/( 2 v/ 2 ), 
while for 77 ^ 1 + k the leading-order behavior is he = rj. The 
important conclusion is that the critical field can be enhanced 
by increasing DMI, thus allowing for faster DW propagation. 


FIG. 6 . (Color online) Critical applied field he vs. DMI parameter 77 . 


IV. DISCUSSION AND CONCLUSIONS 


In recent years there have been ongoing efforts t o use fe r- 
romagnetic materials with perpendicular anisotropjEUllEo] 
produce sharp and stable domain walls with a view to make 
potential spintronic logic and memory devices more com¬ 
pact and faster. Here we have described an alternative ap¬ 
proach to this goal via domain walls in thin nanotubes with 
Dzyaloshinskii-Moriya interaction. These DWs are found 
to have novel properties: The domains themselves become 
twisted about the nanotube, forcing the DWs to become 
sharper with increasing DMI, the opposite of what is seen in 
thin nanowires.ES 

Under applied currents in the regime specified by Eqs. GB 
and ( [T^ , these DWs propagate without distortion with a ve¬ 
locity proportional to the current. Applying a magnetic field, 
we find a rich dependence of DW velocity on polarity, chi¬ 
rality, DMI and field strength, which may provide enhanced 
control in future spintronic devices. The DW velocity can be 
significantly increased by DMI, and the onset of the magnonic 
regime suppressed. 

This work provides the favorable material trends for engi¬ 
neering nanotubes with DMI for faster and more robust DW 
operation. Using the DMI parameter from Ref. [T7l(Z3M? = 
0.5 • 10 “^ J/m^), A = 10 “^^ J/m, and taking the nanotube 
radius R « 100 nm, we estimate the dimensionless DMI pa¬ 
rameter 77 « 2. Eor the same material parameters k = 1 is 
reached for K = 10^ J/m^. These estimates show that the 
regime where DMI has visible effects is experimentally feasi¬ 
ble. 

In the thin nanotube limit, we are able to treat leading con¬ 
tributions of dipolar interactions exactly. We have derived ex¬ 
plicit analytic expressions for the DW profiles and their veloc¬ 
ities under applied currents and fields that are in good agree¬ 
ment with numerical solutions of the LEG equation. These 
results are robust and potentially applicable even beyond the 
thin-nanotube limit, which is hinted at by recent micromag- 
netic studies.® 
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Appendix A: Domain walls of different chirality and polarity 

There are four distinct DW profiles, characterized by po¬ 
larity a — ±1 and chirality x = see Eig.[7] The polarity 
determines whether the DW is head-to-head (cr = 1) or tail-to- 
tail (cr = — 1 ), while chirality determines the sense of rotation 
of m with increasing C, so that sgn 0 q = ax- 


Appendix B: Derivation of Eq. (13) for current driven 
domain-wall motion 

We start with rewriting the LEG equation, Eq. in the 
Landau-Lifshitz (EL) form. Taking the vector product of m 
with both sides of the LEG equation, we obtain 

m X —— = — 7 m X (m X H) 

dt ’ ^ ’ 


dm 9m ^ 

— -Jm X - pj 


dt 


dz 


dz 


(Bl) 


Then, combining LEG equation and Eq. •ED to eliminate the 
m X dm/dt term, we find 

(1 -f — ~ X H — Q!7m X (m x H) 

, , 9m , , 9m 

— (1 -I- al3)J— -(a — j3)Jra x 


dz ’ 
(B2) 


which leads to the EL equation: 
9m 


dt 


= — 7 m X H — am X (m x H) 

~9m - - 9m 

-f J— -f /3Jm X — , 
dz dz 


where 


7 = 


7 


1 -I- a^ 


a = ap = 

J = -' + 
P = 


ap 


l + a 2 ’ 
aj 3 


1 -f a2 
a — d 
1 + a/3 


J, 


(B3) 

(B4) 

(B5) 

(B 6 ) 

(B7) 


^=^/ = (l + a^)^/. 
27 A ^ ' 27 A 


(B 8 ) 


the EL equation takes the form 


9m , / , N 9m 9m 

— =-mxh-amx(mxh)-fji—-fj 2 mx —, (B9) 
dr 9C 9C 


where 


MgR J ! ^MgR 
71 = + 


2A 7 
MgR /3 J 


27 A 


(BIO) 


MgR 

andh= ( = z/R. 

Taking e = “C 1 and ji,j 2 — 0(1), we have that 

the radial component of m is 0(e), and m lies nearly tangent 
to the nanotube. Let us decompose h into its tangential and 
normal components, h = ht -f hn, where ht = h — (h • ep)ep 
and hn = (h • ep)ep. The EL equation (B9 1 is now equivalent 
to a system of two equations: 


9m , / . N dm 

^—=-m X hn - am X (m X ht)-f ji—- , (B12) 

dr dQ 


0 = —m X ht — am x (m x hn) + x 


9m 


(B13) 


Equation (|B12) is the projection of Eq. (|B9) on the tangent 
space of the cylinder, and Eq. ( |B13 1 is the projection of 
Eq. (B91 on Op direction. Erom Eq. (B13 1 we obtain 


, . . 9m 

ahn = m X ht - 32m x —— 

9C 


and consequently. 


u 1 / u ^ , 72 9m 

m X hn = —m X (m x ht) H-^ . 

a a dC, 


(B14) 


(B15) 


The substitution of Eq. ( |B15 1 into Eq. ( |B12) yields 

j 2 \ 9m 


9m / 1 ', 

_ = _ U + - ) m X (m X ht) 


= - a -f 


7i 


a J d( 


m X (m X ht) + j 


9m 


where 


J = - a -f - 
a 


72 \ J 


Rescaling time once again. 


r = I a H— 
a 


-1 


(B16) 


(B17) 


(B18) 
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FIG. 7. (Color online) Domain wall profiles in a thin nanotube with Dzyaloshinskii-Moriya interaction and polarity a and chirality x- The 
magnetization lies tangent to the nanotube surface. 


so that 


may write Eq. equivalently as 


t 


= a 


MsR^ 

2xA 


(B19) 


we obtain 


9m , , , 9m 

— =-mx(mxht)-j —. 


(B20) 


This is the central equation, Eq. ( [T3| , that we analyze through¬ 
out the rest of the paper. 


Appendix C: Expansion of domain wall velocity in magnetic 
field 


+ V-heCOS'd) = {j-v)p. (Cl) 

In terms of the mechanical analogy of Eig. this corresponds 
to energy balance. Letting e = we expand p = Po + ePi + 
e^P 2 and u = uq -F eui + e^V 2 - At zeroth order we deduce 
thatpo(t?(C)) = 0o(C)j where 9q is the static profile, i.e. given 
by Eq. Q. Then it follows that uq = j. The equations for 
the next two corrections pi and p 2 can be readily solved, and 
vi and V 2 are then obtained by integrating Eq. ( |C1| ) over the 
interval 5 < t? < tt + (5 and noting that po vanishes at the 
endpoints. Up to terms of order /ig/a®, we obtain 


While Eq. ( [T8| cannot be solved analytically, it is straight- + k + a'^ tti] 

forward to develop an expansion in powers of hg. It turns out v = j + a -F ^e- iC2) 

that quadratic order is sufficient to capture the leading depen¬ 
dence on polarity and chirality. Letting p(i9(()) = we 
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